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Using the ground-state energy of 16 obtained with the realistic Vucom interaction as a test case, we present 
a comprehensive comparison of different configuration interaction (CI) and coupled-cluster (CC) methods, an- 
alyzing the intrinsic advantages and limitations of each of the approaches. In particular, we use the importance- 
truncated (IT) CI and no-core shell model (NCSM) schemes with up to 4-particle-4-hole (4p4h) excitations, 
with and without the Davidson extensivity corrections, as well as the size extensive CC methods with a com- 
plete treatment of one- and two-body clusters (CCSD) and a non-iterative treatment of connected three-body 
clusters via the completely renormalized correction to the CCSD energy defining the CR-CC(2,3) approach, 
which are all capable of handling larger systems with dozens of explicitly correlated fermions. We discuss 
the impact of the center-of-mass contaminations, the choice of the single-particle basis, and size-extensivity 
on the resulting energies. When the IT-CI and IT-NCSM methods include the 4p4h excitations and when the 
CC calculations include the lplh, 2p2h, and 3p3h clusters, as in the CR-CC(2,3) approach, we observe an ex- 
cellent agreement among the different methodologies, particularly when the Davidson extensivity corrections 
are added to the IT-CI energies and the effects of the connected three-body clusters are accounted for in the 
CC calculations. This shows that despite their individual limitations, the IT-CI, IT-NCSM, and CC methods 
can provide precise and consistent ab initio nuclear structure predictions. Furthermore, the IT-CI, IT-NCSM, 
and CC ground-state energy values obtained for 16 are in reasonable agreement with the experimental value, 
providing further evidence that the Vucom two-body interaction may allow for a good description of binding 
energies for heavier nuclei and that all of the methods used in this study account for most of the relevant particle 
correlation effects. 



I. INTRODUCTION 

Recent developments in nuclear structure theory have lead 
to new perspectives for a QCD-based ab initio description of 
p-shell nuclei and beyond. Significant advances have been 
made regarding the two main ingredients for an ab initio treat- 
ment of nuclear structure: the realistic interaction defining the 
relevant Hamiltonian and the solution of the quantum many- 
body problem. In addition to traditional realistic interactions, 
including the Argonne VI 8 U] or the CD Bonnpotential J2j], 
the framework of chiral effective field theory yl m nas been 
used to construct consistent two- and many-nucleon interac- 
tions based on the relevant degrees of freedom and the sym- 
metries ofQCD HSU. 

Based on these bare interactions, several methods have 
been developed to adapt the nuclear Hamiltonian to the lim- 
ited model spaces available in practical many-body calcu- 
lations. Different conceptual frameworks ranging from the 
Vi ovi k renormalization group method 11811 and the similarity 
renormalization group scheme | ^|^|Tlll to the unitary corre- 
lation operator method lfl2l,[l3l 14lll5ll are being used to con- 
struct soft, phase-shift equivalent interactions, which exhibit 
superior convergence properties when employed in many- 
body calculations. 

Bare or transformed interactions are the starting point for 
modern computational approaches for the ab initio solution of 



the nuclear many-body problem. Two methods have become 
particularly successful in this area, namely, the Green's func- 
tion Monte Carlo approach lfl^ll^[l8lll9l 2fjl and the no-core 
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BoTEl. For nuclei up to A w 12, 
detailed studies of ground and low-lying excited states, in- 
cluding detailed spectroscopic information, were conducted 
with both methods using different two-nucleon interactions 
and phenomenological three-nucleon forces. Very recently, 
the first ab initio no-core shell model calculations using the 
chiral two- plus three-nucleon interactions have shown the 
grea t potential that this new category of interactions may offer 

Although the Green's function Monte Carlo and no-core 
shell model approaches are successful in describing nuclei 
up to the mid-p-shell region, very few methods can provide 
a computationally tractable ab initio description of nuclei 
such as le O or beyond. Among those are the importance- 
truncated no-core shell model and the related importance- 
truncated configuration interaction methods, and the coupled- 
cluster approach, which have recently been used for nuclear 
structure calculations in the A = 16 and A = 40 mass regions 
IHSBIHSSSEE^^SIIIIEI- Extension of 
these methods toward heavier nuclei is presently in progress, 
within both the ab initio description and the traditional effec- 
tive Hamiltonian approach (cf. Refs. 11531 154TJ ). 

The coupled-cluster and configuration interaction method- 
ologies that have been used in these recent developments are 
complementary regarding their respective advantages and lim- 
itations. This work aims at a direct and quantitative compari- 
son of the results of large scale ab initio importance-truncated 
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configuration interaction, importance-truncated no-core shell 
model, and coupled-cluster calculations based on the same 
realistic Hamiltonian using the example of the le O ground 
state. We focus on the most practical coupled-cluster and con- 
figuration interaction approaches that have been specifically 
designed to handle larger systems with dozens of correlated 
particles and larger single-particle basis sets. The most es- 
sential formal differences, advantages, and limitations of the 
importance-truncated configuration interaction, importance- 
truncated no-core shell model, and coupled-cluster methods 
are addressed and their relevance for practical nuclear struc- 
ture calculations is discussed. The nucleon-nucleon interac- 
tion used in this work is the Vucom potential derived in the 
Unitary Correlation Operator Method (UCOM) d H H 
Q 

Following a description of the importance-truncated con- 
figuration interaction and coupled-cluster methods in Sees. [TT] 
and|IIIl we discuss the benchmark results for le O in Sec. HV] 
We assess in detail the impact of center-of-mass contamina- 
tions and the role of the single-particle basis, and compare the 
results of large-scale calculations based on the importance- 
truncated configuration interaction, importance-truncated no- 
core shell model, and coupled-cluster methodologies with one 
another and with experiment. 



II. IMPORTANCE-TRUNCATED CONFIGURATION 
INTERACTION 

As a first class of methods, we consider different configu- 
ration interaction (CI) approaches. The common element of 
all CI methods, which are frequently used in nuclear structure 
theory and quantum chemistry, is that the eigenvalue problem 
for the Hamiltonian is solved in a many-body basis of Slater 
determinants or symmetry-adapted configuration state func- 
tions constructed from a given single-particle basis. In par- 
ticular, the so-called full CI approach of quantum chemistry 
employs a model space spanned by all possible Slater deter- 
minants built from a finite set of single-particle orbits, pro- 
viding the exact solution of the Schrodinger equation in that 
single-particle basis. Unfortunately, the dimension of the full 
CI model space grows factorially with the numbers of parti- 
cles and single-particle orbits, limiting full CI calculations to 
small systems and relatively small single-particle basis sets. 

Truncated CI models offer a simple and natural way of 
reducing the dimension of the model space and the pro- 
hibitive costs of full CI calculations. In the simplest scheme 
that defines the single-reference truncated CI models, only 
the excited determinants up to and including m-particle-m- 
hole (mpmh) excitations from the reference determinant \ $o), 
where m is typically much smaller than the number of active 
particles, are included in the Hamiltonian diagonalization. By 
introducing the fcpfch excitation operators Ck, the ground state 
emerges then as 

m 

|*o> = ^C fe |$ ) , (1) 

fc=0 



where Co $0) is the reference determinant contribution to 
|\I/o) and the Ck |$o) terms with k — 1, . . . , m are the fcpfch 
components of the wave function included in the calculations. 
The popular approaches in this category of methods include 
CISD (CI singles and doubles), CISDT (CI singles, doubles, 
and triples), and CISDTQ (CI singles, doubles, triples, and 
quadruples), if we use the quantum-chemistry acronyms, or 
Cl(mpmh) with m = 2, 3, and 4, respectively, if we use 
the particle-hole language of many-body physics. Since these 
truncated CI methods lack size extensivity (i.e. unlinked di- 
agrams are present in the truncated CI wave function expan- 
sions, resulting in a potential loss of accuracy as the system 
becomes larger), the Davidson extensivity corrections 15511 are 
usually added to the resulting energies to alleviate the prob- 
lem. 

Another popular scheme to define truncated CI models, 
which is frequently used in quantum chemistry and which is 
related to the importance-truncated nuclear CI methods con- 
sidered in this work, is to consider the lplh and 2p2h excita- 
tions from the multi-configurational reference space spanned 
by a number of determinants that provide a reasonable zero- 
order description of the many-body state or states of inter- 
est, as in the multi -reference CI methods lf56l l57ll . The multi- 
reference CI approaches are particularly useful when a single 
reference determinant is a poor approximation to the eigen- 
state(s) of interest and when one wants to accelerate conver- 
gence toward the results of the exact diagonalization by se- 
lecting the dominant contributions to the wave function due to 
3p3h, 4p4h, and other higher-order excitations. With the ju- 
dicious choice of the multi-determinantal reference state, the 
multi-reference CI methods can be very accurate, but they of- 
ten lead to long wave function expansions. Thus, to further 
reduce computer costs of multi-reference CI calculations, one 
often considers reduced sets or subsets of the resulting many- 
body basis states, which are usually defined through the inter- 
nal contractions of configuration state functions or the multi- 
configurational perturbation theory analysis coupled with nu- 
merical thresholds to reject unimportant configurations. The 
importance-truncation scheme used in this work uses similar 
ideas to those exploited in the multi-reference CI approaches, 
while enabling one to incorporate the higher-order mpmh ex- 
citations from the multi-determinantal reference space in a 
systematic fashion. Again, since the multi-reference CI meth- 
ods lack size extensivity, the multi-reference generalizations 
of the Davidson extensivity correction J57l |58l |6(1 |6ll 
are usually added to the resulting energies (for the alternative, 
more intrinsic ways of approximately restoring size extensiv- 
ity in multi-reference CI calculations, see, e.g., Ref. ll62ll and 
references therein). 

In nuclear structure theory, the CI concept is widely used 
in the form of the diagonalization shell model and all of its 
variants (see Ref. ll4lll and references therein). We only con- 
sider no-core calculations in this work, i.e. all particles are 
active. A full CI calculation would use a complete model 
space spanned by all Slater determinants that one can ob- 
tain for a given finite set of single-particle states, e.g. the 
harmonic-oscillator single-particle states up to a maximum 
principal quantum number e max , where e = 2n + I, with n 
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and I representing the radial and angular momentum quan- 
tum numbers, respectively. Other finite sets of single-particle 
states, e.g. those extracted from a self-consistent Hartree-Fock 
calculation, can be employed in the CI calculations as well. 
The nuclear truncated Cl(mpmh) approaches, where m < A, 
emerge then by considering the suitable subsets of Slater de- 
terminants from the full CI model space, including the refer- 
ence determinant |$o) and all kpkh excitations from |$o) with 
k = 1, . . . , m, as described above. By systematically increas- 
ing m, one approaches the full CI limit in which m = A, and 
by considering the e max — > oo limit, one approaches the exact 
solution of the Schrodinger equation for a given Hamiltonian. 
In the following, we will reserve the term CI exclusively for 
the diagonalization shell model calculations based on this spe- 
cific form of the single-particle truncation. 

A modification of the above full and truncated CI meth- 
ods, which has been used in several highly successful ab 
initio investigations of nuclei in recent years, is the no-core 
shell model JNCSM) f~~ 
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formalism uses a different definition of the complete model 
space than that employed in full CI. Instead of considering 
all Slater determinants that one can obtain for a given set 
of the harmonic-oscillator single-particle states up to a maxi- 
mum principal quantum number e max , in the NCSM approach 
one includes those Slater determinants with an unperturbed 
excitation energy that does not exceed N max hil, where is 
the freq uency associated with the harmonic oscillator basis 

Mj^UMjgijM m m mi us m m in m m 



^^EBUi^liliTThe exact solution of the Schrodinger 
equation for a given Hamiltonian is systematically approached 
by considering the N max — » oo limit. Although the NCSM 
calculations rely on the general CI concept of the diagonaliza- 
tion of the Hamiltonian used in all shell-model calculations, in 
this work we distinguish between methods that employ all or 
some Slater determinants resulting from the truncation of the 
single-particle space, which we continue calling the CI meth- 
ods, and the NCSM approaches that represent the special type 
of shell-model calculations which uses the N max hn model 
spaces. Typically, for the reasons explained below, the NCSM 
calculations are performed in the harmonic-oscillator bases 
(cf., however, Ref. for the Hartree-Fock based NCSM 
study). 

For the application to nuclei, the NCSM model spaces offer 
a few important advantages over the model spaces used in CI 
calculations. Unlike in quantum chemistry, where the issue of 
center-of-mass does not exist due to the Born-Oppenheimer 
approximation, the proper description of self-bound nuclei re- 
quires that special care is taken regarding the center-of-mass 
contamination of intrinsic states. Only the NCSM model 
spaces, in which the cutoff for the Slater determinants in- 
cluded in the calculations is defined by setting up the max- 
imum unperturbed excitation energy at N max M}, allow for 
an exact separation of intrinsic and center-of-mass motions, 
which is necessary to obtain translationally invariant intrinsic 
states, as long as a harmonic-oscillator single-particle basis is 
used. The CI truncations, in which the cutoff for the Slater de- 
terminants included in the calculations is typically defined by 



the number of major harmonic oscillator shells corresponding 
to the maximum principal quantum number e max , as defined 
above, violate the separation of intrinsic and center-of-mass 
motions, resulting in the center-of-mass contaminated states. 
The degree of center-of-mass contamination goes down as 
emax ~~ * 00 an d 171 —> A, but one cannot eliminate it math- 
ematically, as in the NCSM case. Furthermore, the truncation 
with respect to the unperturbed excitation energy used in the 
NCSM calculations automatically identifies the most relevant 
Slater determinants in the wave function expansion, as is clear 
from elementary perturbative arguments. One of the questions 
addressed in this work is to examine the potential impact of 
the center-of-mass contamination in truncated CI calculations 
on the results for le O, particularly when compared to other 
sources of errors that all approximate many-body methods us- 
ing computationally tractable truncated model spaces carry. 

For practical applications in nuclear structure calculations, 
the dimensions of the NCSM and CI model spaces become 
prohibitively large, even in relatively small systems, including 
the le O nucleus examined in this work, and even when we use 
the truncated CI approaches. For example, in the aforemen- 
tioned CISDTQ = CI(4p4h) method one has to deal with 
huge numbers of ~ n^n^ and ~ n^n^ determinants of the 
3p3h and 4p4h types and the iterative diagonalization steps 
that scale as n\rft u , where n and n u are the numbers of 
occupied and unoccupied single-particle states, respectively. 
Therefore, in this study we introduce an additional physi- 
cally motivated truncation of the CI and NCSM model spaces, 
which defines the importance-truncated (IT) CI and NCSM 
methods |H1 |H. The IT-CI and IT-NCSM schemes are 
similar in the overall philosophy to the configuration selec- 
tion techniques used in the quantum-chemical MRD-CI model 
of Peyerimhoff and Buenker J57l l64i l65ll . Using the multi- 
configurational perturbation theory of Ref. ll66ll . we estimate 
the importance of individual basis states for the description of 
a certain target state. The amplitude with which a many-body 
basis state contributes in the first-order perturbation the- 
ory to an initial approximation I^Vef) of the target state of 
interest, e.g. the ground state \^o), defines an a priori impor- 
tance measure 



(2) 



where e„ — e re f is the unperturbed excitation energy result- 
ing from a M0ller-Plesset-type partitioning of the Hamiltonian 
116611 . With the intrinsic Hamiltonian H mt = T — T cm + Kjcom 
being a two-body operator, the non-zero importance measures 
can be obtained only for the determinants |$„) that differ 
from the reference state l^ref) by a 2p2h excitation at most. 

In order to access higher-order mpmh excitations with 
to > 2, we embed the concept of the importance measure into 
an iterative construction of the importance-truncated model 
space. Starting from the lowest-energy Slater determinant 
|$o) as an initial reference |^ re t), all basis states with the 
importance measure Ik^I > K m i n , where K m j n defines the im- 
portance threshold, are included in the model space and the 
eigenvalue problem is solved. The resulting eigenvector, pro- 
viding an improved approximation for the target state, is used 
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to define the reference state l^ref), which is now a multi- 
determinantal state, for the second iteration and the above 
procedure is repeated. In the limit « m i n — * 0, this iterative 
scheme converges to the full CI or NCSM space. The pertur- 
bative character of the importance measure entails an mpmh 
hierarchy. In the first iteration, only the determinants up to 
the 2p2h level are generated, resulting in the IT-CI(2p2h) ap- 
proximation. In the limit K m ; n — > 0, the IT-CI(2p2h) calcula- 
tion recovers the full CISD = CI(2p2h) model space. In the 
second iteration, up to 4p4h determinants are present, leading 
to the CISDTQ = CI(4p4h) model space in the n min — > 
limit. In this study, we restrict ourselves to two iterations of 
IT-CI and IT-NCSM, i.e. model spaces up to the 4p4h level, 
so that the resulting IT-CI(4p4h) wave function expansions are 
approximations to the CISDTQ wave functions. We should 
emphasize that for le O and for larger single-particle basis sets 
used in this study, the full CISDTQ = CI(4p4h) calculations 
would not be possible due to the prohibitively large dimen- 
sions of the corresponding model spaces. The IT methodol- 
ogy is among very few methods that enable one to incorporate 
4p4h excitations in the CI calculations for nuclei with larger 
numbers of correlated particles, such as le O, and large single- 
particle basis sets with hundreds of orbitals in a computation- 
ally tractable fashion. 

Eventually, as in the MRD-CI approach of Refs. HI 
l65ll . we construct the importance-truncated model spaces and 
solve the corresponding eigenvalue problems for a sequence 
of importance thresholds /t m j n and extrapolate to K m i n — > 0. 
All IT-NCSM and IT-CI results presented in Sec. [IV]are based 
on sequences of calculations with the importance thresholds in 
the range K min = 2 x 10~ 5 to 15 x 10~ 5 . In order to war- 
rant a robust extrapolation result, we use additional informa- 
tion from a second-order perturbative estimate for the energy 
contributions A e xci(frmin) of the excluded determinants 
i.e. those determinants |$, y ) for which |k„ < K m i n - In the 
limit of a vanishing threshold, which means that no determi- 
nants are discarded, we obviously have A exc i(0) = 0. We 
employ this property as a constraint in a simultaneous extrap- 
olation of the energies E(K m i n ) obtained from the truncated 
eigenvalue problem and of the perturbatively corrected ener- 
gies E(n m i n ) + A exc i(«; m i n ) using 5th to 7th order polynomi- 
als for each of the quantities. The details of this procedure 
will be discussed elsewhere l63ll . The uncertainties of the ex- 
trapolated energies reported in Sec. [IV] are always below 0.5 
MeV. 

The process of constructing the importance-truncated 
model spaces, combined with the extrapolation to the K m ; n — > 
limit, as described above, is very efficient in reducing the 
dimensions of the corresponding CI or NCSM spaces to the 
many-body basis states that are most relevant for the eigen- 
value problem of interest. Nevertheless, we still have to 
solve relatively large eigenvalue problems. Since the resulting 
eigenvalue problems involve sparse Hamiltonian matrices, we 
can use the Lanczos or Arnoldi technique to solve them. We 
use the ARPACK library for that purpose. The largest dimen- 
sions of the importance-truncated model spaces that appear 
in the calculations for le O presented in Sec. [IV] are a factor 
of a few times 10 7 . Since we solve the eigenvalue problem 



in a restricted space, the importance-truncated CI and NCSM 
calculations fulfill the variational principle and the Hylleraas- 
Undheim theorem lloTll . Furthermore, as in all CI and NCSM 
calculations, the IT-CI and IT-NCSM approaches provide us 
with easy access to the eigenstates of the Hamiltonian in a 
shell-model representation, along with the energies, without 
any additional effort. The IT-CI and IT-NCSM eigenstates 
can readily be used to compute expectation values of various 
observables, transition matrix elements, or densities and form- 
factors. 

If we stop the iterative construction of the importance- 
truncated model space before full self-consistency is reached, 
e.g. after just two iterations, as is done in this work, then 
the resulting IT-CI(4p4h) approach is not strictly size exten- 
sive. Thus, as in multi-reference CI calculations in quantum 
chemistry, where a similar problem occurs, we employ the 
aforementioned multi-reference Davidson (MRD) correction 
in order to estimate the effect of higher-order configurations 
beyond the 4p4h excitation level and to approximately restore 
size extensivity. We use the multi -reference rather than single- 
reference Davidson correction, since, in analogy to the multi- 
reference CI methods of quantum chemistry, the IT-CI(4p4h) 
wave function expansions contain subsets of the 3p3h and 
4p4h determinants resulting from the up to 2p2h excitations 
from the multi-configurational reference state ^Vet'}- There 
are several ways of calculating the MRD extensivity correc- 
tions 03, US m US HH IH. We use the Davidson-Silver 
form of the MRD correction discussed in Ref. 16811 . which 
we add to the IT-CI(4p4h) energy to obtain the final, approx- 
imately size extensive IT-CI(4p4h)+MRD result and which is 
calculated using the energy of the multi-configurational ref- 
erence state l^ref) and the summed weight of the reference 
determinants in the IT-CI(4p4h) eigenstate. By employing the 
a posteriori MRD corrections to IT-CI(4p4h) energies, we can 
incorporate the effects of the most essential higher-than-4p4h 
configurations in a computationally efficient manner without 
dealing with the higher-than-4p4h excitations explicitly. In 
the coupled-cluster methods discussed in Sec. [Ill] such con- 
tributions are represented by the disconnected product terms 
involving one- and two-body clusters. As shown in this work, 
the approximately size extensive IT-CI(4p4h)+MRD approach 
is competitive with the rigorously size extensive and accurate 
coupled-cluster schemes that are discussed in the following 
section. 



III. COUPLED-CLUSTER METHOD 

The second class of methods that we consider are those 
based on the single-reference coupled-cluster (CC ) theory 
M M HD, H! H (see Refs. [zl M WlM M 
for selected reviews), which utilizes the exponential ansatz for 
the A-particle ground state, 

|*o) = exp(T)|$ ), 0) 
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where |$o) is the reference determinant and 

A 

fe=i 

is the cluster operator. The cluster operator T is a particle- 
hole excitation operator, defined relative to the Fermi vacuum 
| $0)1 whose many -body components 

Tk= £ %;:Xal---ala ik --- ail (5) 

il<---<ik,a 1 <---<a k 

generate the connected wave function diagrams of |^o) to in- 
finite order. The remaining linked, but disconnected contri- 
butions are produced through the exponential ansatz for l^o)- 
Here and elsewhere in this paper, we use the notation in which 
ii, - ■ ■ or i, j, . . . are the single-particle states occupied in 
the Fermi vacuum state |$ ) an d a i, a 2, ■ ■ • or a, b, . . . are the 
single-particle states unoccupied in |$o) (we will use labels 
p,q, . . . for the generic single-particle basis states). 

Typically, the explicit equations for the ground-state energy 
Eq and the cluster amplitudes Vg'" 1 * defining the many -body 
components Tk of T according to Eq. (0 are obtained by in- 
serting the wave function l^o), Eq. (0, into the Schrodinger 
equation, H\^o) = -Eb |^ r o)» premultiplying both sides of the 
resulting equation on the left by cxp(— T) to obtain the con- 
nected cluster form of the Schrodinger equation l7lH72il . 

H\$ ) =Eo\*o), (6) 

where 

H = exp(-T)ifexp(T) = [H cxp(T)] c (7) 

is the similarity-transformed Hamiltonian or, equivalently, the 
connected product of the Hamiltonian and exp(T) (desig- 
nated by subscript C), and projecting both sides of Eq. (0 
on the reference determinant |$o) an d excited determinants 
= < ■ • ■ ■ ■ • OijSo) that span the relevant 

y4-particle Hilbert space. The latter projections result in a non- 
linear system of explicitly connected and energy-independent 
equations for the cluster amplitudes t^"*£ , 

|# |$ ) = 0, h < ■ ■ ■ < i k , 01 < • • • < o fc> (8) 

where H is defined by Eq. (0 and k = 1, . . . , A, whereas the 
projection on |$o) results in the CC energy formula, 

E Q = ($ |#|$o>. (9) 

The advantage of this formulation of CC theory, which is a 
standard formulation adopted, for example, in quantum chem- 
istry, where most of the development and application work 
involving CC methods has occured, is that, unlike the ex- 
pectation value of the Hamiltonian with the CC wave func- 
tion, which would lead to a non-terminating power series in 
T of the form E = ($ \[exp(T*)H exp(T)] c |$ ) JH, 
the CC equations written above represent algebraic expres- 
sions that mathematically terminate at a finite power of T. 



The similarity-transformed Hamiltonian H, or the connected 
product of the Hamiltonian and the exp(T) operator which is 
equivalent to H [cf. Eq. (0], is a finite polynomial expansion 
in T whose length depends only on the highest many-body 
rank of the interactions in the Hamiltonian, not on the num- 
ber of particles in a system, so one does not need to make 
ad hoc truncations in powers of T, which would have to be 
invoked if one attempted to minimize the expectation value 
of the Hamiltonian, in order to solve for the cluster ampli- 
tudes and energy. For example, for Hamiltonians that do not 
contain higher-than-pairwise interactions, H and the result- 
ing CC equations for cluster amplitudes, Eq. (0, terminate at 
the T 4 terms, which is a purely mathematical truncation re- 
sulting from the fact that one cannot connect more than four 
T diagrams to the diagrams representing the Hamiltonian in 
the definition of H , Eq. dT), if H is a two-body Hamiltonian. 
The energy formula, Eq. (0, simplifies even further in this 
case. If H does not contain higher-than-pairwise interactions, 
we obtain 

E Q = ($ |i/|$o) + ($o|[ J ffw(Ti+T 1 2 + iT 1 2 )] c |$ ), (10) 

where iZjv = H — (<&q\H\$q) is the Hamiltonian in the 
normal-ordered form relative to |$o)- In other words, after 
determining the cluster operator T by solving the system of 
equations represented by Eq. ([8]), we only need the lplh or 
singly excited and 2p2h or doubly excited components of T, 
T\ and Ti, respectively, to determine the energy if H is a two- 
body Hamiltonian. 

The above is the exact CC theory, which is equivalent to the 
exact diagonalization of the Hamiltonian with the full CI ap- 
proach. Indeed, we could obtain Eqs. (© and (0 by directly 
projecting the Schrodinger equation for the exact CC wave 
function |* ), Eq. ©, Hexp(T)|$ ) = E exp(T)|$ ), 
with T defined by Eq. ©, on the |$ > and l^ 1 //.^} de- 
terminants that span the A-particle Hilbert space. The energy- 
dependent terms on the right-hand side of the resulting equa- 
tions, 

(S>%;;%\HeMT)\*o) =^o(^ 1 1 ;.t fc |exp(T)|$ > ! (11) 

cancel out the unlinked terms on the left-hand side of Eq. ( ITTb 
to produce the system of the explicitly connected amplitude 
equations represented by Eq. (0 (see, e.g., Ref. 08111 for a 
pedagogical derivation). Similarly, the disconnected diagrams 
corresponding to the product of the Hamiltonian and cxp(T) 
in the resulting energy formula, 

E a = <$o|#exp(T)|$ }, (12) 

do not contribute to Eq. ( fT2l . resulting in Eq. (0. Alterna- 
tively, one can show that Eq. (0 is equivalent to the expec- 
tation value of the Hamiltonian with the CC wave function, 
($o|exp(rt) jffe xp(r)|$o}/($o|exp(rt)exp(T)|$ ), as 
long as the cluster operator T has the exact form given by 
Eq. (0, in which all many-body components of T including 
Ta are included, and as long as the corresponding cluster 
amplitude s ffi '"'*, with k = 1, . . . , A satisfy Eq. (0 (see, 
e.g., Ref. Mi" " 
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The exact CC theory, as described above, being equivalent 
to the full CI diagonalization, is limited to small few-body 
problems. Thus, in all practical applications of CC theory, 
including those reported in this work, one truncates the many- 
body expansion for T at some, preferably low, mpmh excita- 
tion level T m . In this study, following the footsteps of quan- 
tum chemistry, where the CC theory has become one of the 
most successful and frequently used many-body methodolo- 
gies, we focus on the most practical CC approximations that 
can be applied to systems with dozens or even hundreds of 
correlated fermions. Thus, we consider the basic CCSD (CC 
singles and doubles) approximation ll83l [84l l8~5l l86ll . which 
accounts for the effects of one- and two-body clusters, T\ 
and T 2 , respectively, as well as the completely renormalized 
CR-CC(2,3) approach JH H IH, which accounts for the 
effects of connected three-body T3 clusters through a rela- 
tively inexpensive, yet very effective non-iterative correction 
to the CCSD energy and which represents an improved vari- 
ant of the cornpletelyrenormalized CCSD(T) [CR-CCSD(T)] 
method lf77ip78i[8lll90ll used in previous ab initio no core cal- 
culations for the 4 He and le O nuclei SSHHHS. In addi- 
tion to many successf ul q uantum chemistry applicat ions (see 

Refs. niisiiiiSiiiSSmiiaEilili^ifor 

representative examples), the CR-CC(2,3) approach was re- 
cently used to study the ground and excited states of the 56 Ni 
nucleus, treated by the effective Hamiltonian in the p/-shell 
basis l53ll . demonstrating its ability to provide a virtually ex- 
act description of closed-shell nuclei at a tiny fraction of the 
costs of the shell-model calculations that aim at similar ac- 
curacies, such as the CI method with up to 4p4h excitations, 
and the ability to retain the accuracy of the CI approach with 
up to 4p4h excitations even when the reference determinant 
contributes very little to the wave function [the CPU timings 
of the CR-CC(2,3) calculations reported in Ref. ll53tl were on 
the order of one minute on a single processor]. The present 
paper shows the high accuracy offered by the CR-CC(2,3) ap- 
proach in the context of the much larger scale no-core ab ini- 
tio calculations for le O employing realistic nucleon-nucleon 
interactions and single-particle bases as large as 8 major os- 
cillator shells (480 uncoupled single-particle basis states). 

In the CCSD calculations, which, in addition to produc- 
ing reasonably accurate results for closed-shell systems, pro- 
vide the framework for the determination of the CR-CC(2,3) 
corrections due to T3 clusters, the many-body expansion of 
the cluster operator T defining the CC wave function ansatz 
is truncated at the level of 2p2h (or double) excitations, i.e. 
T = T 1 + T 2 , where [cf. Eq. ©] 

Ti=^Ao,i (13) 

i . a 

and 

r 2 = £ C44%^ (W) 

i<j,a<b 

are the singly and doubly excited clusters. In analogy to the 
exact CC theory defined by Eqs. (0, where k = 1, . . . , A, and 
©, we determine the cluster amplitudes tf a and defining 



the CCSD wave function 

|*( CCSD) }=exp(T 1+ T 2 )|$ } (15) 

by solving the nonlinear system of energy-independent and 
explicitly connected algebraic equations similar to Eq. (0. 
Specifically, since we only need the Ti and T 2 cluster com- 
ponents, we consider the subset of equations corresponding to 
the projections on the singly and doubly excited determinants 
[Eq. ([8} with k = 1 and 2], so that the number of equations 
matches the number of unknown amplitudes t l a and t^ b . We 
obtain 

($«|£-(ccsd)| $ ^ = Qj (16) 
($«6|jj(CCSD)| $( ^ =Q) 

where 

g(ccsv) = cxp (_r 1 -T 2 )i/exp(r 1 +T 2 ) 

= [frexp(Ti+T 2 )] (18) 

is the similarity-transformed Hamiltonian of CCSD [the 
similarity-transformed Hamiltonian H, Eq. (|7]), written for 
T = Ti + T 2 ], and and |$° h ) are the singly and doubly 
excited determinants, respectively, relative to |$o)- Once the 
cluster amplitudes t % a and t** b defining T\ and T 2 are known, 
the CCSD energy Eq CSB ^ is computed using Eq. (O in 
which H is replaced by ( CCSD ), Eq. jTSJ. 

As pointed out above, the CC equations obtained through 
projections of the Schrodinger equation have finite polyno- 
mial nature, i.e. they terminate at a finite power of T which 
depends on the highest many-body rank of the interactions in 
the Hamiltonian. When the Hamiltonian is a two-body oper- 
ator, the CC equations for cluster amplitudes terminate at the 
T 4 terms and the ground-state energy is given by Eq. ( TTOb . In 
particular, when H does not contain higher-than-two-body in- 
teractions, the generic CCSD amplitude equations, Eqs. ( fT6l ) 
and ( fTTI i, that apply to all Hamiltonians simplify to 

($?|[lf J v(l+T 1 +T 2 + |T 1 2 +T 1 T 2 + iT 1 3 )] c |$o) = (19) 

and 

{$$\[H N (1 +Tx + T 2 + \Tl + TiT 2 + \Tf + ±T| 

+^i 2 T 2 + ^T 1 4 )] c |$o) = 0, (20) 

respectively. If the cluster operator T was not truncated, as 
in the exact CC theory defined by Eqs. ©, ©, and ©, the 
results would be equivalent to the variational, full CI calcu- 
lation corresponding to the minimization of the expectation 
value of the Hamiltonian with the CC wave function |^o)> Eq. 
(fJJ. Unfortunately, in the CCSD case, T is truncated at T 2 , so 

( CCSD) 

that the resulting ground-state energy Eq , determined 
by solving the CCSD amplitude equations, Eqs. (TToT l and ( TT7t 
[or, when the Hamiltonian is two-body, (fT~9b and d20l)l, and 
by using the resulting Ti and T 2 clusters in Eq. (O [or, in 
the two-body Hamiltonian case, (TTOb l is not equivalent to the 
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<5o(2,3) 



E 



i j k 
be 



(21) 



i<j<k,a<b<c 



to the CCSD energy E { Q CCSD \ so that the final CR-CC(2,3) 



£0t 



ibe 



expectation value of the Hamiltonian with the CCSD wave energy is 

(CCSD) 

function }, Eq. (1151 1. and, as such, does not have the 

bound if A > 2. One could, at least in principle, contemplate 
variational CCSD calculations based on minimizing the ex- 
pectation value o f the Ham iltonian with the CCSD wave func- 
tion, as in Refs. |l0lUl02ll . but, as already explained, the ex- 
pectation value of the Hamiltonian with the CC (e.g. CCSD) 
wave function is a non-terminating series in cluster amplitudes 
that does not lead to practical computational schemes. On 
the other hand, due to the exponential nature of the CC wave 
function, the approximate CC methods using truncated cluster 
operators T = J2T=i Tk with m < A (CCSD corresponds to 
the to = 2 case) converge more rapidly to the full CI limit as 
m — ► A than the equivalent truncated shell-model expansions 
that rely on the same manifold of excited determinants, de- 
fined by Eq. (Q]), while providing the rigorously size extensive 
results which truncated CI expansions cannot provide due to 
the presence of unlinked contributions in the wave function 
and disconnected contributions in the energy that do not can- 
cel out. For example, as already established in the 1970s and 
early 1980s, the CCSD method, which relies on up to 2p2h 
cluster components, is more accurate than the CI approach 
truncated at 2p2h determinants. As mentioned earlier, and as 
demonstrated in this article, the CCSD approach corrected for 
the connected 3p3h clusters T3 via the CR-CC(2,3) method 
is as accurate as the CI approach with up to 4p4h excitations, 
while providing a size extensive description of many-particle 
systems. The issue of the lack of extensivity in CI calcula- 
tions is addressed in this paper through the suitably defined, 
quantum-chemistry inspired corrections to the CI energies, as 
described in the previous section. The elementary analysis 
that enables one to understand why the CR-CC(2,3) method 
is as accurate as the CI approach with up to 4p4h excitations 
is provided in Section JVC] 

The CCSD approximation provides information about the 
bulk of the correlation effects, which are described by T\, T 2 , 
and their products, but one must account for the higher-than- 
doubly excited connected clusters, particularly the T3 clusters, 
in order to obtain a more quantitative description. This can be 
done through the CCSDT (CC singles, doubles, and triples) 
method 11031 110411 . in which the cluster operator is truncated 
at the 3p3h level (T = T x + T 2 + T 3 ) and, in addition to 
the equations corresponding to the projections on singly and 
doubly excited determinants, one considers the projections on 
triply excited determinants [k = 3 in Eq. dSJ]. Unfortu- 
nately, as explained below, the CCSDT scheme has very large 
computational costs, and so is generally not practical and lim- 
ited to small few-body problems and relatively small single- 
particle basis sets. As a result, inspired by the similar diffi- 
culties encountered in quantum chemistry, where CCSDT is 
limited to small few-electron systems, we instead use the con- 
siderably less expensive CR-CC(2,3) approach, in which we 
incorporate the effects of T3 clusters by adding a non-iterative 
a posteriori correction 



E, 



(CR-CC(2,3)) 




E 



(CCSD) 




<*o(2,3), 



(22) 



where So(2, 3) is given by Eq. (fJTJ. The quantities 9JtS, en- 
tering Eq. (fJTJ are the projections of the connected cluster 
form of the Schrodinger equation written for the CCSD wave 
function on the triply excited determinants which 



-ijk I 

define t he triply exci t ed moments of the CCSD equations 

i.e. 



define the triply excited mo 

fnl EzEIHl Ha, M M 0, l 



(23) 



where _ff( CCSD ) is the similarity-transformed Hamiltonian of 
the CCSD approach, Eq. (Q~8]l. As in the case of the CCSD 
equations, the explicit form of Eq. ( l23l depends on the nature 
of the interactions in the Hamiltonian. When the Hamiltonian 
does not contain higher-than-two-body interactions (the case 
examined in this paper), one can write 



S$g\[H N (T2 + T X T 2 + ±T 2 2 + |Tf3 

■\t x tI + irfr 2 )] c |$ ). 



The coefficients entering Eq. (f2TTi are defined as 



nabc 



(<f> Q \AHW%)/D*< 



abc\ 1 j^ijk 
be' 



(24) 



(25) 



where the denominator 



D 



ijk 
abc 



? (CCSD) 
v 



$ o6c| 5 (CCSD)| $ «6c > (26) 



ijk 1 



is obtained by approximating the triples-triples block of the 
matrix representing the i//( CCSD ) operator, Eq. dT8l . by its 
diagonal, as in the Epstein-Nesbet perturbation theory, while 
A = Ai + A 2 is the hole-particle deexcit ation operator defin- 
ing the CCSD 'bra' or dual ground state ifiol [TM [l07ll 



(CCSD) I 




($o|(l + Ai+A a )exp(-Ti-r 2 ), (27) 



which satisfies the biorthonormality condition 



(CCSD)|^(CCSD) 











} = 1. The one- and two-body compo- 



nents of the A operator of CCSD, 



and 



A 2 = K] a \ a ] a baa, 

i<j,a<b 



(28) 



(29) 



respective ly, are obtaine d by solving the linear system of 
equations 

(<P \{l + A 1 +A 2 )H^ cs ^\^) =4 CCSD) A?, (30) 



($o|(l + Ai + A 2 )i?( CCSD ) |*#) = E { 



(CCSD) , Qh 







(31) 
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The details of the derivation of the non-iterative correction 
5o(2,3) defining the CR-CC(2,3) calculations, which orig- 
inates from the mathematical formalism referred to as the 
biorthogonal method of moments of coupled-cluster equations 
J87i[88ll and which describes the leading terms toward the full 
CI energy due to T3 clusters, mimicking the performance of 
the much more expensive CCSDT method at the small frac- 
tion of the computer cost, can be found in Refs. I87ll88ll . 

The fact that the construction of the CR-CC(2,3) correction 
($o(2,3) requires the determination of the A operator of the 
CCSD theory, which defines the CCSD bra state (^ CCSD) |, 
as described above, has an additional advantage that we can 
use the same operator A to determine properties other than 
energy. We can, for example, use it to determine the CCSD 
one-body reduced density matrices 

Tj = (§ CCSD) |(qt a9 )|^r SD) ) 

= ($o|(l + Ai+A 2 )(4a,)|$o) (32) 

and properties 

(#( CCSD )|9|* GCSD) ) = £0*7?, 03) 

where is a one-body property operator defined through ma- 
trix elements 9^ = (p\9\q), 



(apa g ) = exp(-Ti - r 2 )(aj,a g )exp(Ti + T 2 ) 

= [{ala q )cMTi+T 2 )] c (34) 

is the similarity-transformed connected form of the operator 
string (ajpdq), analogous to the similarity-transformed Hamil- 
tonian are the singly and doubly excited 

clusters obtained by solving the CCSD equations, Eqs. (TToT l 
and (T% , and Ai and A 2 are obtained by solvin g the linear sys - 
tem of equations given by Eqs. (f30| and ED IU0 l[l06l[l07ll . 
In general, if 5 is an operator representing the quantum- 
mechanical quantity of interest, the CC analog of the ex- 
pectation value of 3 can be determined using the equation 

(5) = (*o|S|* ) = (*o|(l + A) 5|$ ), (35) 

where 

3 = exp(-T) 3 cxp(T) = [H cxp(T)] c (36) 

is the similarity-transformed form of 3 and A is the hole- 
particle deexcitation operator defining the left ground eigen- 
state of H, Eq. ©, which has the form ($ (1 + A), and the 
bra or dual CC ground state 

(*o| = (*o|(l + A)eaq)(-r). (37) 

In analogy to the CCSD case, A can be obtained by solving 
the linear system of equations, 

($ |(i + A)#|$£v;£ } = p: \?; :: % , (38) 



where Eq is the CC ground-state energy and A? 1 "'"* are the 
amplitudes that define the many-body components of A, 

Afe= E X tt k 4 1 ---4 k ^ k ---a ai (39) 

il<---<jfc,ai<---<a fc 

[clearly, Eq. d38l ) is the generalization of Eqs. d30l ) and ED* 
corresponding to the CCSD case, to any level of CC theory]. 
In particular, when 3 is the Hamiltonian H , we obtain from 
Eq. E3 

(H) = <$ |(l+A)iT|$o>, (40) 

which is equivalent to the CC energy formula, Eq. (O, since 
the (& \A k H |$ ) contributions to Eq. d40b vanish when the 
cluster operator T entering the definition of H , as in Eq. (0, 
satisfies the system of CC equations, Eq. ([8]). It should be 
noted that the above way of calculating the (3) values, which 
reflects on the biorthogonal character of CC theory, becomes 
fully equivalent to the determination of (3) as the conven- 
tional expectation value (5 , o|3| v I / o}/(^'o| v I'o} when T is a 
non-truncated cluster operator given by Eq. obtained in 
the exact CC calculations defined by Eq. (E)). When T is 
truncated, as in the CCSD case, the value of (3) determined 
from Eq. E3, although no longer identical to the quantum- 
mechanical expectation value (\I> |S|\I/o) / (^0 1 ^0 ) ^ where 
l^o) is the corresponding CC wave function, is equivalent to 
the alternative way of determining (5) as (dEo(X)/d\)\ = o, 
where Eg(X) is the CC energy, Eq. (O, calculated after solv- 
ing the relevant CC equations for the auxiliary Ha milto nian 
H\ = H + A3, as in the response CC theory JH Il08ll . as 
long as the reference determinant |$o) does not vary with A. 
We use this fact in Section llV Al to determine the expectation 
values of the center-of-mass Hamiltonian corresponding to the 
CCSD and CR-CC(2,3) calculations as the appropriate CC en- 
ergy derivatives. The vast experience with performing CC cal- 
culations in quantum chemistry is that the difference between 
the values of (3) calculated as the conventional quantum- 
mechanical expectation values and the (3) values obtained as 
the corresponding energy derivatives, as described above, are 
very small, since the approximate CC methods, such as those 
used in this work, provide results close to full CI. Clearly, 
there are no differences between the (3) values obtained as the 
traditional expectation values and energy derivatives in trun- 
cated CI calculations, since all CI methods are variational and, 
as such, satisfy the Hellmann-Feynman theorem. In approx- 
imate CC methods, we have to rely on a response formula- 
tion and the equations such as Eq. d35l l. or the equivalent en- 
ergy derivatives, as described above, since, in analogy to the 
Hamiltonian, the traditional expectation value expression with 
the CC wave function, (\E , |3|* )/( v I , o|^ , o), would lead to a 
non-terminating power series in cluster amplitudes of the form 
(i> | [exp(T t ) 3 exp(T)]c| < l ) o} that does not lead to practical 
computational schemes. 

By using the CR-CC(2,3) approach, we are able to signifi- 
cantly reduce the computational costs associated with the in- 
clusion of the connected triply excited clusters, enabling cal- 
culations with relatively large model spaces, while producing 
results that should be as accurate as those obtained with the 
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full CCSDT scheme (§7l]. Indeed, the most expensive compu- 
tational steps of CR-CC(2,3) scale as n^nf t in the determina- 
tion of the non-iterative correction due to T3 and n^n^ in the 
underlying CCSD calculation, where, as mentioned earlier, 
n a and n u represent the numbers of occupied and unoccupied 
single-particle states, respectively. For realistic values of n 
and n u , including the larger single-particle basis sets used in 
this work, this is less expensive than the n^n^ iterative steps 
defining CCSDT by orders of magnitude. In addition, unlike 
in CCSDT, in the CR-CC(2,3) calculations one does not have 
to store the large number of ~ n^n 3 amplitudes t l L defin- 
ing the T3 cluster operator, since the non-iterative correction 
Jo(2, 3) defining the T3 correction to the CCSD energy, Eq. 
(f2TT >. is calculated using the one- and two-body clusters, T\ 
and T2, and their A deexcitation analogs, Ai and A2, respec- 
tively, as described above, which need a storage of the n n u 
amplitudes t l a and Xf and n^rv^ amplitudes t^ b and Xfj. Thus, 
the total storage requirements of the CR-CC(2,3) calculations 
are primarily defined by the number of two-body matrix ele- 
ments that define the Hamiltonian or the two-body matrix el- 
ements of the similarity-transformed Hamiltonian H ( CCSD ) . 
The memory requirements scale as n n 3 . A small subset of 
the three-body matrix elements of H^ CCSI) > that enter the cal- 
culation of the CR-CC(2,3) triples correction Jo (2, 3) via the 
aforementioned equations for the A operator of CCSD and 
the triply excited moments •EDT^ 7 fec do not have to be precom- 
puted and stored, since one can rigorously factorize the dia- 
grams that represent them and express all quantities that enter 
the calculation of Jo (2, 3) in terms of the o ne an d two-body 
matrix elements of fl"( CCSD ) (see, e.g., Ref. [ 109], and refer- 
ences therein). Thanks to the consistent use of the recursively 
generated intermediates, including one- and two-body matrix 
elem ents of ( CCSD ), and fast matrix multiplication routines 
ill 1 Oh - as is always done in the most efficient, modern imple- 
mentations of CC methods in quantum chemistry, our CCSD 
and CR-CC(2,3) com puter codes are fully vectorized. Thanks 
to the DIIS algorithm ill ill , which we use to solve the CCSD 
equations and their A counterparts (for the first application 
of th e DIIS procedure to solving the CC equations, see Ref. 
ill 12|P . our CCSD calculations typically converge in about 20 
iterations to obtain an energy accurate to within 10~ 5 MeV. 
The triples correction Jo(2, 3) requires no iterations, which is 
one of the biggest advantages of CR-CC(2,3) as opposed to 
the iterative CCSDT approach. We refer the reader to Refs. 
Hpj [H for the details of the CR-CC(2,3) theory, Ref. 
109h for the explicit, computationally efficient expressions 
for the one- and two-body matrix elements of H ( CCSD ) and 
moments ^JCL that enter the CR-CC(2,3) calculations, and 



Ref. H5 lH for the computationally efficient, fully factorized 
form of the CCSD equations that precede the calculation of 
the triples correction Jo(2, 3). We also note that while the 
earlier form of the CR-CC(2,3) approach, designated as CR- 
CCSD(T) ||^|^[8l)j9pLwhich was used in the studies of the 



*He and lb O nuclei [|47| 



, 48t], was only approximately size ex- 
tensive (to within 0.5-1 % of the correlation energy ll77lo . the 
CR-CC(2,3) method used in this work is fully size extensive, 
so that no loss of accuracy occurs when going from smaller to 
larger many-fermion systems. 



One of the most important findings of this work is that the 
CR-CC(2,3) method, with its relatively inexpensive computa- 
tional steps that scale as n^n^ in the iterative part and v? n\ 
in the non-iterative part, is capable of providing the results 
that are virtually identical to those that effectively correspond 
to the diagonalization of the Hamiltonian using the CI ap- 
proach with up to 4p4h excitations from the reference deter- 
minant I $0)1 corrected for the effects of higher-than-4p4h ex- 
citations and extensivity through the use of the Davidson cor- 
rections. The IT-CI(4p4h) approach reduces the most expen- 
sive n^n^ steps of the full CI(4p4h) calculation and the need 
to deal with the large numbers of the 3p3h and 4p4h determi- 
nants in the CI(4p4h) diagonalization by orders of magnitude 
through the use of importance truncation and extrapolation to 
the K m in — + limit, as discussed in the previous section. The 
CR-CC(2,3) method does effectively the same work, eliminat- 
ing, in particular, the need to deal with the 3p3h and 4p4h ex- 
citations of CI in an explicit manner, by representing the dom- 
inant higher-order correlations through the disconnected prod- 
uct terms involving low-order T\ and T2 clusters, such as, for 
example, T\T 2 for 3p3h excitations, (1/2)T 2 2 for 4p4h excita- 
tions, (l/2)TiT 2 2 for 5p5h excitations, (1/6)T 3 for 6p6h exci- 
tations, etc., while making sure that the connected 3p3h exci- 
tations that are usually more important than the disconnected 
T1T2 terms are accounted for through the computationally af- 
fordable non-iterative Jo (2, 3) correction to the CCSD energy. 
The detailed analysis of the CI and CC wave functions that ex- 
plains this is provided in Sec. IIVCI The numerical similarity 
of the CR-CC(2,3) and Davidson-corrected IT-CI(4p4h) [i.e. 
IT-CI(4p4h)+MRD)] results does not address the issue of the 
center-of-mass contaminations which are present in all trun- 
cated CC and all truncated CI calculations other than NCSM. 
The impact of the center-of-mass contaminations on the trun- 
cated CC [CCSD and CR-CC(2,3)] and C I [IT-CI(4p4h)] re- 
sults for le O is examined in Section lTV AI 



IV. BENCHMARK FOR 16 

We aim at a quantitative comparison of the different ap- 
proaches using the ground state of le O as an example. For 
the interaction we use the realistic Vucom potential derived 
in the framework of the Unitary Correlation Operator Method 
(UCOM) discussed in Refs. (Qj, GJ, 0, [H ■ Using a unitary 
transformation tailored for the description of short-range cen- 
tral and tensor correlations, a phase-shift equivalent soft inter- 
action is derived from the Argonne V18 potential |1:]. Appli- 
cations in various many-body approaches, from the no-core 
shell model for light nuclei to Hart ree-F ock and many-body 
perturbation theory for heavy nuclei 111 1311 . show that this two- 
body potential allows for a realistic description of binding en- 
ergies throughout the nuclear chart without the explicit inclu- 
sion of additional three-body (or other higher-than-two-body) 
interactions. All calculations presented here are based on the 
sam e Vu com (I0 = 0.09 fm 3 ) potential as that used in Refs. 

02 El. 

As implied by the discussion in the previous section, we 
compare the following three groups of ab initio methods: 
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(i) the coupled-cluster approach including singly and dou- 
bly excited clusters, CCSD, as well as the non-iterative com- 
pletely renormalized CR-CC(2,3) scheme that corrects the 
CCSD results for the effects of the connected triply excited 
clusters, (ii) the importance-truncated configuration interac- 
tion method including up to 4p4h excitations, IT-CI(4p4h), 
without and with the MRD size extensivity corrections, and 
(iii) the importance-truncated no-core shell model approach 
including up to 4p4h excitations, IT-NCSM(4p4h), with and 
without the MRD corrections. In addition to comparing the 
results of CC, CI, and NCSM calculations, particularly the 
CR-CC(2,3) and IT-CI(4p4h)(+MRD) levels, we examine sev- 
eral issues relevant for this comparison, including the role of 
center-of-mass contaminations and the sensitivity of the re- 
sults to the choice of the single-particle basis. By comparing 
the size extensive, but not variationally bound CR-CC(2,3) re- 
sults with the variational, but not rigorously size extensive IT- 
CI(4p4h)(+MRD) results for the binding energies of le O side 
by side, we have an opportunity to comment on the signifi- 
cance or insignificance of such issues as the violation of size 
extensivity by the truncated CI calculations and the lack of the 
variational bound in the CC calculations. 

We begin with two important issues relevant for the com- 
parison of the CCSD, CR-CC(2,3), IT-CI(4p4h)(+MRD), and 
IT-NCSM(4p4h) results, namely, the role of center of mass 
contaminations and the sensitivity of the results to the choice 
of the single-particle basis. 



A. Center-of-mass problem 

All calculations are performed with a translationally invari- 
ant intrinsic Hamiltonian H mt = T — T cm + Vucom- However, 
this does not imply that the resulting many-body states and 
intrinsic energies are free of spurious contributions induced 
by a coupling of intrinsic and center-of-mass (CM) degrees 
of freedom. For the nucleus as a self-bound and translation- 
ally invariant system the intrinsic properties should not de- 
pend on the CM motion, i.e. intrinsic and CM components of 
the many-body state have to decouple. 

For a Slater determinant basis, an exact separation of 
intrinsic and CM motions is possible only in a complete 
]V max /jO model space based on harmonic oscillator single- 
particle states as employed in the full NCSM approach. The 
use of a different model-space truncation or a different single- 
particle basis destroys this decoupling property and induces 
CM contaminations. This problem is well k nown in the con- 
text of the nuclear shell model (see Ref. Il 14fl and refer- 
ences therein) and it was also addressed in the context of 
CC calculations for nuclei, both intrinsical ly, through a trans- 
lationally invariant formulation 11 15l 11 ldl . and numerically, 
through heuristic CM corrections added to the Hamiltonian 

JH SI SI m H m, Hi in m . 

In order to probe to what extent intrinsic and CM motions 
are coupled in the different many-body approaches examined 
in this work, we study the impact of an artifi cial shift of th e 
CM spectrum using the Lawson prescription 11 14lll 1 7L 1 1 18ll . 



TABLE I: The CI and NCSM expectation values of H^t and H cm , and 
their CC analogs defined in the text (in units of MeV), obtained from 
the many-body solutions using Hp for (3 = and (3 = 10. The CC 
and IT-CI calculations use a harmonic oscillator basis with e max = 5 
and HQ = 30 MeV, whereas the IT-NCSM and NCSM calculations 
use a model space with iV max = 8 and hQ, = 22 MeV. The oscillator 
frequencies correspond to the respective energy minima of the CC 
and NCSM calculations. 
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with CM momentum P cm and CM coordinate X cm . If the CM 
motion is completely decoupled, then the expectation value 
of the intrinsic Hamiltonian H mt , (H mt ), computed with the 
ground state resulting from the modified Hamiltonian Hp is 
independent of the value of (3. Any dependence of (iJ; n t) on (3 
indicates an unphysical coupling of intrinsic and CM motions 
and a violation of translational invariance. 

In Table |T] we compare the results for the binding energy 
of le O obtained with the different many-body methods used 
in this work for two (3 values, namely, (3 = 0, at which the 
initial intrinsic Hamiltonian is recovered, and j3 = 10, which 
is a t ypical empirical value used in shell model applications 
Il 14ll . We use the representative Ml values which approxi- 
mately match the minima on the curves that describe the de- 
pendence of the relevant energies on Ml Since there is no rig- 
orous criterion for choosing (3, we will come back to the im- 
pact of variations of this parameter later on. The IT-CI(4p4h), 
IT-NCSM(4p4h), and NCSM values of (H int ) and (H cm ) are 
calculated directly as the standard quantum-mechanical ex- 
pectation values of the relevant Hamiltonians using the eigen- 
vectors obtained from the solution of the eigenvalue prob- 
lem of Hp. The CCSD and CR-CC(2,3) analogs of (H mt ) 
at (3 = are calculated using the appropriate CC energy for- 
mulas [Eq. OH for CCSD and Eq. <22) for CR-CC(2,3)] 
applied to H = H mt , rather than the expectation values of 
H mt with the CC wave functions which, as explained in Sec- 
tion [nfl are not used in the practical implementations of CC 
methods employed in this work. The CCSD and CR-CC(2,3) 
(H cm ) values are obtained by numerically differentiating the 
corresponding energies, Eq. ( fTOb for CCSD and Eq. ( 1221 for 
CR-CC(2,3), where H — Hp, with respect to (3 at the (3 values 
of interest, as in the response CC considerations described in 
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FIG. 1: (color online) The intrinsic ground-state energies, (-Hint), of 
16 0, obtained with the CCSD and CR-CC(2,3) approaches at fi = 
(open symbols) and f3 = 10 (filled symbols), as functions of ftfi for 
a harmonic oscillator basis with e max = 5: CCSD (•.o) and CR- 
CC(2,3) (♦,<». 



the previous section. Following the Lawson recipe, the CCSD 
and CR-CC(2,3) values of (flint) at nonzero (3 are calculated 
as E (/3) - /3(H cm )i3, where E (/3) is the relevant CC [CCSD 
or CR-CC(2,3)] energy obtained for the CM-corrected Hamil- 
tonian Hp and (H cm )p = 3Eq((3) / 3(3 is the corresponding 
value of (H cm ) calculated at the same j3 through the numeri- 
cal differentiation of Eo(f3), as described above. 

As Table U shows, the full NCSM approach allows an ex- 
act separation of intrinsic and CM motions and consequently 
the intrinsic energy (flint) is independent of (3. The IT- 
NCSM(4p4h) method shows a minimal variation of the in- 
trinsic energy at the level of 0.1 MeV, which is a conse- 
quence of the importance truncation. This tiny coupling in 
the IT-NCSM(4p4h) calculations can safely be neglected. The 
CCSD, CR-CC(2,3), and IT-CI(4p4h) methods violate trans- 
lational invariance from the beginning through the choice of 
the model space. This results in a more sizable coupling be- 
tween intrinsic and CM motions when compared to the IT- 
NCSM(4p4h) scheme: the intrinsic energies (flint) change by 
2 to 3 MeV when going from (3 = to f3 = 10. As men- 
tioned above, the oscillator frequencies used in Table U cor- 
respond to the approximate positions of the minima on the 
curves that show the dependencies of the relevant energies on 
Ml. The full foil-dependence of the intrinsic energies (flint) 
for 13 = and = 10 obtained in the CCSD and CR-CC(2,3) 
calculations is depicted in Fig. Q] Evidently, the change of 
(flim), when going from f3 = to f3 = 10, increases some- 
what with increasing oscillator frequency. For example, for 
hfl = 18 MeV, the intrinsic energy obtained with the CR- 
CC(2,3) approach changes by 1.18 MeV when going from 
P = to (3 = 10, whereas the analogous change in the CR- 
CC(2,3) energy for Ml = 34 MeV is 2.03 MeV. 

The above results demonstrate that the Lawson prescription 
is quite efficient in reducing the expectation value of the CM 
Hamiltonian, (fl cm ), obtained from the CCSD, CR-CC(2,3), 
and IT-CI(4p4h) solutions, as shown in Table J] However, this 
does not mean that a decoupling between intrinsic and CM 



motions or a spuriousness-free intrinsic state has been com- 
pletely achieved. In fact, if we take the same Lawson prescrip- 
tion and compare the intrinsic energies (flint) for f3 = 10 and 
(3 = 20, we observe changes in the values of (flint) which are 
similar in magnitude to those observed when one goes from 
f3 = to (3 = 10. For example, for the IT-CI(4p4h) calcu- 
lations at hfl = 30 MeV using 6 major oscillator shells (i.e. 
emax = 5)> m e change in the (flint) value when going from 
(3 = to 13 = 10 is 1.35 MeV. When going from /3 = 10 to 
[3 = 20, the IT-CI(4p4h) value of the (fl int ) energy changes 
by 0.91 MeV although (fl cm ) at {3 = 10 is already below 0.2 
MeV. This shows that the smallness of (fl cm ) alone does not 
indicate a decoupling or warrant an intrinsic state free of CM 
contaminations. 

Apart from an explicit projection | 114|| or a translationally 
invariant formulation, there is no rigorous way to eliminate 
the CM contamination problem from the CC and CI calcu- 
lations. Out of the methods presented in this work, only the 
NCSM and IT-NCSM calculations are free (in the IT-NCSM 
case, virtually free) of this spuriousness. The CC and CI re- 
sults are CM contaminated although, as shown in Tableland 
Fig. [T] the degree of this contamination, when the harmonic 
oscillator reference |$o) is used, seems relatively small in the 
case of le O. Indeed, the degree of CM contamination, as mea- 
sured by the changes in the intrinsic energies (fl n t) when go- 
ing from (3 = to f3 = 10, does not seem to exceed 2-3 
MeV when a medium-size single-particle basis consisting of 
6 major oscillator shells is employed and |$o) is the harmonic 
oscillator reference, although we must remember that simi- 
lar changes in the intrinsic energies are observed when go- 
ing from f3 = 10 to (3 = 20, suggesting that the real degree 
of CM contamination in the CC and CI results is somewhat 
bigger than 2-3 MeV. The degree of CM contamination, as 
measured by the (fl cm ) values at (3 = 0, does not seem to ex- 
ceed about 5-6 MeV in the same basis. When the size of the 
single-particle basis is increased, as is done in the following 
sections where we examine single-particle basis sets as large 
as e max = 7 (8 major oscillator shells), the magnitude of the 
unphysical coupling of intrinsic and CM motions in the CC 
and CI calculations is expected to be reduced, since we ap- 
proach the limit of the complete single-particle basis. On the 
other hand, other effects, such as the ^-dependencies of the 
resulting energies shown in Fig. Q~]and the fact that the CCSD, 
CR-CC(2,3), and IT-CI(4p4h) methods do not provide the ex- 
act wave function that would factorize into the intrinsic and 
translational components, might hinder the reduction of the 
CM contaminations present in the CC and CI results. For all 
these various reasons, in the assessment of the quality of the 
CC and CI calculations reported in this work, we will remain 
cautious and keep in mind that the resulting ground-state ener- 
gies may carry an uncertainty anywhere between 2 and 6 MeV 
or so as a result of CM contamination, at least when the har- 
monic oscillator reference |$o) is employed. We will continue 
examining the role of CM contaminations on the CCSD, CR- 
CC(2,3), and IT-CI(4p4h) calculations with different types of 
single-particle bases and different basis set sizes in the future 
work. 
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B. Role of the single-particle basis 

In many cases, the harmonic oscillator basis is not the opti- 
mal choice for the expansion of the nuclear many-body state. 
The naive reference determinant |$o} obtained by occupying 
the lowest-energy harmonic oscillator states may have a rel- 
atively small overlap with the final wave function re- 
sulting in unnecessarily long CI expansions to represent the 
correlated |^o) state that have to compensate for the defi- 
ciencies of the reference determinant | <&o ) - By switching to a 
single-particle basis optimized for the nucleus under consider- 
ation, generated, for example, by a Hartree-Fock calculation, 
the convergence with respect to the many-particle basis can 
be significantly enhanced. Generally, this option is not used in 
the NCSM calculations, since the use of a single-particle basis 
set other than the harmonic oscillator basis would destroy the 
mathematical decoupling of intrinsic and CM motions that the 
NCSM model space guarantees for any finite basis set. How- 
ever, for the truncated CC or CI approaches, where we do not 
have this property anyway, we may benefit from the use of 
optimized single-particle bases. 

In order to demonstrate the effect of an optimization of the 
single-particle basis, we compare the results of the CC and CI 
calculations using the harmonic oscillator (HO) and Hartree- 
Fock (HF) bases. The latter is obtained from a self-consistent 
Hartree-Fock calculation using the same intrinsic Hamiltonian 
-Him and the same single-particle space as tho se employed in 
the subsequent many -body calculations ill 13ll . Therefore, the 
HF optimization can be viewed as a unitary transformation 
within the set of single-particle states employed in the many- 
body calculations. For a full CI calculation at given e max > 
where all determinants resulting from a given single-particle 
basis set are included, this transformation would not affect the 
results. The situation changes when one uses the truncated 
CI and CC wave function expansions, where the HO — > HF 
transformation of single-particle states may have an effect on 
the resulting energies. As elaborated on below, this effect is 
expected to be small in the case of CC calculations, which rely 
on the exponential form of the wave function that makes the 
results of truncated CC calculations virtually invariant with re- 
spect to orbital rotations through the presence of the exp(Ti) 
com ponen t in the CC wave operator, as in the Thouless theo- 
rem ill 19ll . but can be quite significant when the truncated CI 
expansions are employed, since the CI wave operator is a lin- 
ear rather than an exponential excitation operator which does 
not contain sufficiently many terms to make the results numer- 
ically invariant with respect to orbital rotations if truncated at 
the mpmh excitations with m <C A. Our numerical analysis 
confirms these expectations. 

The impact of the basis optimization on a IT-CI(4p4h) cal- 
culation is illustrated in Fig. [2 a), where we compare the HO- 
and HF-based IT-CI(4p4h) results obtained with the e max = 5 
model space. We performed similar calculations using other 
model spaces up to and including e max = 7, and the results 
are very similar to those shown in Fig.|2a), so in the following 
discussion we focus on the e max = 5 case. At smaller oscil- 
lator frequencies, the ground-state energies of le O obtained 
with both bases agree very well. However, with increasing 
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FIG. 2: (color online) Ground-state energy of 16 O as a function of hQ, 
obtained with the HO single-particle basis (open symbols) and the 
HF-optimized basis (filled symbols) in an e max = 5 model space, 
(a) Importance-truncated configuration interaction calculations: IT- 
CI(4p4h) («,0) and IT-CI(4p4h)+MRD (♦,<». (b) Coupled-cluster 
calculations: CCSD («,0) and CR-CC(2,3) (♦,<»• 



Kl the HF basis leads to lower ground-state energies than the 
HO basis. At the same time, the Davidson correction for the 
HO-based calculation increases rapidly, indicating that contri- 
butions beyond the 4p4h level become significant in this case. 
In contrast to the HO-based IT-CI(4p4h) calculation, the HF- 
based IT-CI(4p4h) calculation develops a minimum at larger 
Kl and the Davidson correction remains small at all frequen- 
cies, clearly implying that the role of higher-than-4p4h excita- 
tions in the CI expansion of the ground-state wave function of 
le O is suppressed when the optimum HF determinant is used 
as a reference determinant | <&o) - 

The analogous analysis for the CCSD and CR-CC(2,3) cal- 
culations, presented in Fig. |2jb), reveals that the CC meth- 
ods are significantly less sensitive to the choice of the single- 
particle basis. Again, we only show the sample of the HO- and 
HF-based CCSD and CR-CC(2,3) calculations correspond- 
ing to the e max = 5 model space. We performed similar 
CC calculations using other model spaces up to and includ- 
ing e max = 7, and the observed patterns are similar to those 
shown in Fig.|2b), so we focus on the e max = 5 case here. As 
was the case for the IT-CI(4p4h) calculations, the energies ob- 
tained with the two bases are virtually the same for smaller Ml 
values. As we increase the oscillator frequency, the disagree- 
ment between the HO- and HF-based CC results grows, but 
only slightly, particularly for the CR-CC(2,3) approach. In- 
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deed, unlike the IT-CI(4p4h) case, where the largest discrep- 
ancies between the results obtained with the two bases in the 
Ml = 14-34 MeV region, which occur at Ml = 34 MeV, 
are as much as 8 MeV with the Davidson correction and about 
20 MeV without it, the analogous differences between the 
HO- and HF-based CC results at larger Ml values are quite 
small, with the two results differing by 3.8(4.2) MeV in the 
CCSD case and 1.4(1.5) MeV in the CR-CC(2,3) case, when 
the e max = 5 basis set is employed and Ml = 34(38) MeV. 
Furthermore, we see that the overall shapes of the curves dis- 
playing the Ml dependence of the total energies obtained with 
the CCSD and CR-CC(2,3) approaches are very similar re- 
gardless of which single-particle basis is used. It is interest- 
ing to note that for both the CC and CI methodologies, the 
agreement between the results obtained with the two bases 
improves as we use more accurate approximations. This is 
particularly true for the CC calculations, where the differ- 
ence between the HO- and HF-based results obtained with 
the CR-CC(2,3) method, which is a more accurate approx- 
imation when compared to CCSD, is smaller than the analo- 
gous difference between the HO- and HF-based CCSD results. 
The differences between the CI results obtained with the two 
bases are generally larger, but even in this case, the differ- 
ence between the HO- and HF-based results obtained with the 
IT-CI(4p4h)+MRD approach, which corrects the IT-CI(4p4h) 
energies for at least some effects of higher-than-4p4h excita- 
tions, are smaller than the analogous difference between the 
HO- and HF-based IT-CI(4p4h) results. This makes sense, of 
course, since the closer we get to the full CI limit, the less 
sensitive the results become with respect to orbital rotations. 

The relative insensitivity of the CC results to the choice 
of the single-particle basis, which has been kno wn in quan- 
tum chemistry for a long time (cf., e.g., Ref. lll20ll ), is a 
cons equen ce of the implicit inclusion of the Thouless theo- 
rem jl 19ll in the CC calculations through the exp(Ti) com- 
ponent of the CC wave operator exp(T), even when T is 
truncated at the two-body level, as in CCSD, where T = 
T1+T2. The exp(Ti ) component of the CCSD wave operator 
exp(Ti + T 2 ) = exp(Ti)cxp(T 2 ) (Ti and T 2 are particle- 
hole excitation operators and hence they commute), obtained 
by solving the coupled system of equations for the Ti and T 2 
clusters, as described in Section|III] acting on the reference de- 
terminant |$o) effectively optimizes the single-particle basis 
for the many-particle state of interest without the need for the 
explicit introduction of orbital relaxation, thus reducing the 
impact of any inadequacies of the basis on the final results and 
making the CCSD energies almost independent of the type 
of the single-particle basis. One can simply write the CCSD 

wave function, Eq. O, as |<4 CCSD) ) = exp(Ti+T 2 )|$o) = 
exp(T 2 )|$Q), where |$q) = exp(Ti)|$ ) is a new refer- 
ence determinant optimized for the ground-state |'I'q CCSD ' ) ) 
through the suitable choice of T\ obtained from the CCSD 
calculations. Indeed, if we analyze the T\ cluster amplitudes 
resulting from our CCSD calculations, we see that they are 
quite large when the HO basis and large Ml are employed, i.e. 
when the HO reference determinant $0) is far from the op- 
timum reference. When the optimized HF basis is employed, 
the T\ cluster amplitudes resulting from CCSD calculations 



are small, independent of Ml. These patter ns are reflected 
in the values of the so-called 'Ti Diagnostic' 112111 . which is 

defined as ((^oIt/TiI^o}/^) 172 - In defining the T x Diag- 
nostic, we divide the connected (i.e. size extensive) quantity 
($o|T/Ti |$o)i which represents the magnitude of T\ cluster 
contributions to the wave function, by n to make the result 
independent of the system size. For the case of the e max = 5 
model space and the HO basis, the T\ Diagnostic resulting 
from the CCSD calculations changes quite dramatically, from 
0.10 at Ml = 14 MeV to 0.42 at Ml = 38 MeV and 0.48 at 
Ml = 42 MeV, indicating a steep increase in the values of the 
Ti cluster amplitudes due to the increasing inadequacy of the 
HO reference determinant |$o) at larger Ml values that those 
large Ti amplitudes compensate for. For the optimized HF 
basis, the same diagnostic remains almost constant, with val- 
ues that do not exceed 0.04 in the entire Ml = 14 - 42 MeV 
region, confirming the smallness of T\ cluster amplitudes in- 
dependent of Ml in the HF-basis case. In fact, it is easy to un- 
derstand why T\ is small in the HF basis. When | $0 ) is the HF 
reference, the lowest orders of many-body perturbation theory 
that the T\ cluster component contributes to are second for the 
wave function and fourth for the energy. For comparison, T\ 
contributes to the first-order wave function and second-order 
correction to the energy when the non-HF references, such 
as the HO reference determinant, are used. In the region of 
larger Ml values, the HO reference determinant is so far from 
the optimum HF determinant that the resulting T\ amplitudes 
become very large, as reflected in the above values of the T\ 
Diagnostic. 

In order to confirm that the primary role of the T\ operator 
in CC calculations is to effectively relax the orbitals to pro- 
duce the optimum reference determinant |$q) = cxp(Ti)|<l>o) 
for the many-particle state of interest, we compare the T 2 clus- 
ter operators resulting from the HO- and HF-based CCSD cal- 
culations. We expect the T 2 clusters, which describe the lead- 
ing correlation effects, to be very similar in the HO- and HF- 
based CCSD calculations if the main role of T\ is to optimize 
the reference determinant. The 'T 2 Diagnostic' (cf., e.g., Ref. 
JH), which is defined as (($o\T%T 2 \<S>o) /rio) 1 / 2 and which 
measures the significance of the T 2 cluster contributions, con- 
firms this expectation. For the e max = 5 model space, the val- 
ues of the T 2 Diagnostic obtained from the HO-based CCSD 
calculations are 0.17-0.18 in the entire Ml = 14 - 42 MeV 
region. The analogous values of the T 2 Diagnostic resulting 
from the HF-based CCSD calculations are 0.15-0.17 in the 
same region. Thus, the T 2 clusters that describe the true cor- 
relation effects barely change with the type of the basis (HO 
vs. HF) and Ml. The T% clusters remain small and do not 
change much with Ml when the optimized HF basis is em- 
ployed, while becoming sizable in the large Ml region when 
the naive reference determinant |$o)> obtained by filling the 
lowest-energy HO single-particle states, which is a poor rep- 
resentation of the correlated ground state in that region, is em- 
ployed. The CC theory can cope with the inadequacy of the 
HO basis by using the exp(Ti) component of the CC wave 
operator with large Ti amplitudes to transform the naive refer- 
ence determinant $0) resulting from the use of the HO basis, 
as in the Thouless theorem, to the more optimal form which 
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is adjusted to the ground-state wave function \^>o). The same 
exp(Ti) operator does not change the reference determinant 
|$o) too much when T\ is small, i.e. when the orbitals are 
properly optimized beforehand, as in the HF case. This ex- 
plains the virtual invariance of the CC results on the choice of 
the single-particle basis. 

The same arguments enable us to understand why the re- 
sults of truncated CI calculations may significantly depend on 
the type of the basis in the region of larger Ml values and 
why one needs to use the HF-optimized orbitals in that re- 
gion to obtain the results of the CC quality with the trun- 
cated CI approaches. Let us focus on the IT-CI(4p4h) ap- 
proach and the related CISDTQ scheme, in which the linear 
excitation operator C defining the ground-state wave func- 
tion l^o) through the formula l^o) = Cl^o) has the trun- 
cated form C = C + C t + C 2 + C 3 + C 4 . Here, C k 
is the fcpfch excitation operator generating the contributions 
from the fc-tuply excited determinants when acting on |$o) 
(Co generates the reference contribution). When the inter- 
mediate normalization condition (c&ol^o) = 1 is imposed 
on the CI wave function, so that Co becomes a unit opera- 
tor, the lplh component of the CI wave function, Ci\$>o), 
is equivalent to the lplh component of the CC wave func- 
tion, Ti|$o)- In addition to the connected Ti|$o) contribu- 
tion, the IT-CI(4p4h) and CISDTQ wave functions contain 
the disconnected cluster terms, such as (l/2)T 1 2 |<!>o} (through 
the C 2 contribution), (1/6)T 1 3 |$ ) (through the C 3 contribu- 
tion), and (1/24)T 1 4 [$ ) (through the C 4 contribution), but 
they do not contain the entire cxp(Ti)|$o) expansion, which 
includes higher powers of T\ if A > 4 (as is the case for 
the le O nucleus). In other words, the linear excitation opera- 
tors C of the IT-CI(4p4h) and CISDTQ approaches or other 
truncated CI schemes do not have the mathematical struc- 
ture of the Thouless theorem that would enable one to fac- 
tor out the exp(Ti) component which would make the results 
virtually independent of the orbital choice. In consequence, 
the results of truncated CI calculations may display a strong 
dependence on the choice of the basis, as our IT-CI(4p4h) 
calculations shown in Fig. [2] clearly demonstrate. Full CI is 
the only CI method that contains the exp(Ti)|$o) component 
in its entirety, since one can always represent the intermedi- 
ately normalized full CI wave function for the A-body system, 

*o) = (1 + C\ H + Ca)\$o), in the exponential form 

|*o) = exp(T 1 + --- + T A )[$ >. 

The analogous analysis can be used to explain why the 
HO-based IT-CI(4p4h) calculations become less accurate in 
the region of larger Ml values and why the IT-CI(4p4h) cal- 
culations benefit from the use of the HF-optimized orbitals, 
particularly in the region of larger Ml values. As already 
pointed out, the T\ cluster component becomes large in the 
region of larger Ml values when the HO reference |$o) is 
employed. This can be seen by analyzing the CCSD wave 
function, as described above, or by examining the IT-CI(4p4h) 
wave function. Since the lplh component of the CI wave 
function, G\\$>q), is equivalent to the lplh component of the 
CC wave function, Ti|$o)i when the intermediate normal- 
ization is imposed, we can immediately conclude that if T\ 
is large, as is the case for larger Ml values and HO basis, 



the corresponding C\ excitation amplitudes should be large 
as well. This is exactly what we observe in the HO-based 
IT-CI(4p4h) calculations. For example, the largest C\ excita- 
tion amplitude in the IT-CI(4p4h) intermediately normalized 
wave function increases from 0.08 at Ml = 14 MeV to 0.29 
at Ml = 30 MeV and 0.40 at Ml = 38 MeV, when the 
HO basis is used and e max = 5. In contrast, the HF-based 
IT-CI(4p4h) calculations exhibit largest C\ amplitudes which 
do not exceed 0.03 throughout the entire M7-range. Now, if 
the Ci or T\ component is large, the product terms such as 
TxT 2 , (1/2)2?T 2 , (1/6)7^, (l/2)Ti2?, (1/4)2?!?, etc. 
become relatively large as well. All of these terms are in- 
cluded in the CCSD wave function, helping the accuracy of 
the CCSD and CR-CC(2,3) calculations, but not all of them 
are present in the IT-CI(4p4h) wave function, which contains 
the TiT 2 and (1/2)T 1 2 T 2 components through the 3p3h and 
4p4h contributions described by C3 and C4, respectively, but 
not the (1/6)T X 3 T2, (l/2)TiT 2 2 , and (1/4)T 2 T 2 components, 
which represent the 5p5h [(1/6)7^2 and (l/2)TiT 2 2 ] and 
6p6h [(1/4)I?T|] excitations neglected in IT-CI(4p4h). As 
we can see, the absence of the 5p5h, 6p6h, etc. components 
in the IT-CI(4p4h) wave function hurts the accuracy of the IT- 
CI(4p4h) calculations when C\ or T\ is large, which is exactly 
what happens when we use the HO basis in the region of larger 
Ml values. 

The situation dramatically changes when the HF basis is 
employed. In that case, the C\ or T\ contributions are 
small and the higher-order product terms, such as (1/ 6)TfT 2 , 
(l/2)TiT 2 2 , and (1/4)7?T|, which are neglected in the IT- 
CI(4p4h) calculations, become very small as well, resulting in 
an excellent description of the le O nucleus by the IT-CI(4p4h) 
method in the entire region of Ml which matches the accuracy 
of the CR-CC(2,3) calculations, particularly when the already 
very good IT-CI(4p4h) results are corrected for size exten- 
sivity and the remaining small higher-than-4p4h excitations 
through the multi-reference Davidson correction. The excel- 
lent agreement between the HF-based IT-CI(4p4h)+MRD and 
CR-CC(2,3) results, illustrated in Fig. [2] for e max ~ 5, re- 
mains valid when larger model spaces are employed, enabling 
us to draw several important conclusions regarding the qual- 
ity of the ab initio results for le O reported in this work. These 
conclusions are discussed in the next section. 



C. Comparison of large-scale calculations 

We now compare the predictions for the ground-state en- 
ergy of le O obtained in the IT-NCSM, IT-CI, and CC calcula- 
tions employing Vucom and the largest model spaces that we 
could handle with a reasonable computational effort. 

In Fig. |3]the convergence of the IT-NCSM ground-state en- 
ergy as a function of A max is shown for fixed Ml = 22 MeV, 
which is determined by the position of the energy mini- 
mum for the larger model spaces l63ll . In addition to the 
IT-NCSM(4p4h) results without and with the Davidson cor- 
rection, we report the res ults of full NCSM calculations us- 
ing the Antoine code lll22ll for A^ax < 8. In the re- 
gion where these exact reference results are available, the IT- 
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FIG. 3: (color online) Importance-truncated no-core shell model cal- 
culations for the ground-state energy of le O with Vucom for Ml = 
22 MeV. Shown are the results of the IT-NCSM(4p4h) calculations 
(•), IT-NCSM(4p4h)+MRD calculations (♦), and full NCSM calcu- 
lations (+). The vertical line indicates the result of an exponential 
extrapolation of the IT-NCSM(4p4h)+MRD energies. 



NCSM(4p4h)+MRD and full NCSM energies agree to within 
0.1 MeV or better. 

Using exponential fits involving the five largest model 
spaces Il23ll . we obtain an extrapolated ground-state energy 
of (-129 ± 1) MeV for the IT-NCSM(4p4h) data and of 
(-130 ± 1) MeV for the IT-NCSM(4p4h)+MRD results. The 
change due to the Davidson correction provides an estimate of 
the effects beyond 4p4h configurations. This estimate agrees 
very well with the preliminary results of the explicit inclusion 
of up to 6p6h configurations which will be fully elaborated on 
elsewhere 063I1 . The comparison of these re sults with the ex- 
perimental binding energy of —127.6 MeV lll24ll proves that 
the Vucom two-body interaction provides an excellent descrip- 
tion of ground-state energies for heavier nuclei. 

The results of the IT-CI(4p4h) and IT-CI(4p4h)+MRD cal- 
culations based on the HF-optimized basis are summarized 
in Fig. Eta). With increasing size of the single-particle ba- 
sis from e max = 4 to 7 the position of the energy minimum 
shifts systematically towards larger Ml The Davidson cor- 
rection remains on the order of 1 MeV for all model-space 
sizes and oscillator frequencies, indicating that the effects of 
beyond-4p4h configurations are small when the HF basis set 
is employed. A similar picture emerges from the CC calcula- 
tions, shown in Fig. Hfb), which use the same HF-optimized 
bases and the same model spaces as the IT-CI(4p4h) and 
IT-CI(4p4h)+MRD calculations. The inclusion of connected 
triples through the CR-CC(2,3) scheme leads to a lowering of 
the ground-state energy by up to 6 MeV, indicating the impor- 
tance of T3 cluster contributions in the quantitative calcula- 
tions of the nuclear binding energies. 

A direct comparison of the IT-CI(4p4h)+MRD and CR- 
CC(2,3) results is presented in Fig. Hfc). The agreement be- 
tween these two entirely different many-body approaches is 
extraordinary. Apart from the largest model space employed 
in this study consisting of 8 major oscillator shells, the two 
data sets are practically on top of each other. For the largest 
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FIG. 4: (color online) Systematic comparison of IT-CI and CC re- 
sults for the ground-state energy of ie O using HF-optimized single- 
particle bases with e max = 4, 5, 6, and 7. (a) Comparison of IT- 
CI(4p4h) (open symbols) with IT-CI(4p4h)+MRD (filled symbols). 
(b) Comparison of CCSD (open symbols) with CR-CC(2,3) (filled 
symbols), (c) Comparison of IT-CI(4p4h)+MRD (open symbols) 
with CR-CC(2,3) (filled symbols). 



e m ax = 7 space, the discrepancies between the CR-CC(2,3) 
and IT-CI(4p4h)+MRD results are slightly larger than in the 
case of smaller basis sets, i.e. the CR-CC(2,3) energies are 
approximately 1 — 2 MeV lower than the corresponding IT- 
CI(4p4h)+MRD energies, but the overall agreement between 
the IT-CI(4p4h)+MRD and CR-CC(2,3) energies is outstand- 
ing. Based on this systematic agreement, we can conclude that 
neither the lack of strict size extensivity of the truncated IT- 
CI(4p4h) calculations, which can be taken care of through the 
use of the Davidson corrections, nor the violation of the varia- 
tional principle by the CC methods, which is compensated by 
the high accuracy these methods offer, pose significant practi- 
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cal problems. 

The excellent agreement between the CR-CC(2,3) and IT- 
CI(4p4h)+MRD data can be rationalized by comparing the CC 
and CI wave function expansions. If we impose the interme- 
diate normalization condition (c&ol^o) = 1 exploited in CC 
theory on the exact CI wave function expansion, we obtain 
the following relationships between the CI excitation opera- 
tors C n and the CC cluster components T n : 

Ci = Ti, (43) 

C 2 = T 2 + \Tl (44) 

C 3 = Tz+TxT 2 + \Tl, (45) 

d = T 4 + TiT 3 + \Tl + \TlT 2 + ±T?, etc. (46) 

By design, the CR-CC(2,3) approach provides a highly ac- 
curate description of the connected T\, T 2 , and T 3 clusters, 
but not of T4 or T n with n > 4, and of all of the dis- 
connected product terms that enter the C n excitation opera- 
tors with n = 1 - 4, except for T1T3. The T1T3 term is 
much smaller than the leading 4p4h component represented 
by (l/2)Tf, particularly when the HF basis is employed. 
When the HF basis is employed, (l/2)Tf contributes to the 
second-order many -body perturbation theory correction to the 
wave function and the fourth-order correction to the energy, 
whereas the lowest-order corrections to the wave function and 
energy resulting from the T1T3 cluster are fourth and sixth, 
respectively. The connected T4 cluster contributions, which 
contribute to the fifth and higher orders in the many -body per- 
turbation theory expansion for the energy, are much smaller 
than the disconnected (l/2)Xf contributions as well. In fact, 
much of the success of CC theory in areas such as quantum 
chemistry is related to the fact that one can safely neglect T4 
in calculations for the non-degenerate closed-shell systems. 
The negligible role of T4 clusters has also been observed in 
the study of the semi-closed-shell 56 Ni nucleus, as described 
by the effective Hamiltonian in the p/-shell basis l53ll . The 
le O nucleus is a closed-shell system, so one does not need T4 
to accurately describe its ground state. This is why the CR- 
CC(2,3) approach provides a virtually exact description of the 
CI excitation contributions C n up to and including the 4p4h 
(i.e. n = 4) terms. Since the CR-CC(2,3) method is based on 
the idea of correcting the CCSD energy for the leading T3 con- 
tributions and since the CCSD approach describes all higher- 
than-4p4h excitations that can be represented as products of 
the T\ and T 2 clusters, one has to correct the IT-CI(4p4h) en- 
ergies for the selected higher-than-4p4h correlations via the 
Davidson corrections to improve the agreement between the 
IT-CI(4p4h) and CR-CC(2,3) data. This explains why the IT- 
CI(4p4h)+MRD and CR-CC(2,3) results obtained in this work 
agree so well. 

In order to compare the IT-CI(4p4h)+MRD and CR- 
CC(2,3) results with the aforementioned converged IT-NCSM 
calculations, an extrapolation e max — > 00 is necessary. In 
view of the convergence pattern of the CI and CC results and 
the fact that our IT-CI(4p4h)+MRD and CR-CC(2,3) data are 
limited to e max < 7, this extrapolation can only provide a 
rough estimate. The crude exponential extrapolations based 
on the total energies obtained with the four different model 



space sizes with e max = 4 — 7 at fixed ft£l = 30, 34, 38, and 
42MeV, for which the full set of the IT-CI(4p4h)+MRD and 
CR-CC(2,3) data for e max = 4 — 7 is available, lead to an 
estimate of the e max — > 00 CR-CC(2,3) energy in the range 
of -131 to -133 MeV in the entire hfl = 30 - 42 MeV 
region. The IT-CI(4p4h)+MRD result is similar. It looks 
as though the e max — > 00 CR-CC(2,3) energies of —131 to 
-133 MeV and their IT-CI(4p4h)+MRD analogs, obtained via 
the above exponential extrapolations based on the total en- 
ergies, are in excellent agreement with the extrapolated IT- 
NCSM(4p4h)+MRD result of ( - 130±1) MeV, but we should 
be very careful in interpreting this agreement, which might be 
fortuitous, since crude extrapolations based on the limited set 
of CC and CI data that we have access to may carry sizable 
error bars. For example, if we use the more careful approach 
where instead of extrapolating total energies we extrapolate 
only the correlation energies obtained with e max = 3, 5, and 
7, and add the results to the highly accurate estimate of the 
converged (to within 0. 1 MeV) HF energy resulting from the 
e m ax = 20 HF calculations, we obtain the extrapolated to- 
tal CR-CC(2,3) energy in the range of -135 to -141 MeV 
in the entire hfl = 30 — 42 MeV region. The analogous IT- 
CI(4p4h)+MRD energies are in the range of -132 to -135 
MeV. The use of only odd values of e max in the correlation en- 
ergy extrapolations can be justified by the fact that the HF and, 
in consequence, correlation energies do not change uniformly 
when increasing the basis set; changes in the HF energies are 
much stronger when another radial excitation is added to the p 
states, as observed when going from e max = 4 to 5 and from 
emax = 6 to 7. Thus, based on the limited set of the CC and 
CI data we have at our disposal, we can only state that the 
extrapolated e max -> 00 CR-CC(2,3) and IT-CI(4p4h)+MRD 
energies fall within the broader range of —131 to —141 MeV, 
which implies that our extrapolated results carry an uncer- 
tainty which could be as big as about 10 MeV. Clearly, the few 
MeV differences between the extrapolated CR-CC(2,3)/IT- 
CI(4p4h)+MRD and IT-NCSM(4p4h)+MRD results for the 
binding energy of le O prompt further study, so that we can 
understand the nature of these differences in more precise 
terms, but it is already worth pointing out that these few 
MeV differences between the extrapolated CR-CC(2,3) and 
IT-CI(4p4h)+MRD binding energies on the one hand and their 
Nmax — > 00 IT-NCSM(4p4h)+MRD counterpart on the other 
hand are consistent with our estimate of the effect of CM con- 
taminations on the calculated CR-CC(2,3) and IT-CI(4p4h) 
energies discussed in Section HV Al Since the CR-CC(2,3), IT- 
CI(4p4h)+MRD, and IT-NCSM(4p4h)+MRD methods con- 
tain similar many-body correlation effects as e max (7V max ) — > 
00 and since the IT-NCSM(4p4h)+MRD results are virtually 
free of the CM contaminations (cf. Section IIV At . it is quite 
possible that the CM contaminations in the CR-CC(2,3) and 
IT-CI(4p4h)+MRD results are largely responsible for the ob- 
served few MeV differences between the extrapolated CR- 
CC(2,3) or IT-CI(4p4h)+MRD and IT-NCSM(4p4h)+MRD 
energies. The fact that the relatively inexpensive CR-CC(2,3) 
approach can produce the binding energy of le O which dif- 
fers from the best IT-NCSM(4p4h)+MRD estimate and exper- 
iment by only a few MeV (less than 10 %), when the Kjcom 
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interaction is employed, is an indication that the CR-CC(2,3) 
method captures practically all correlations relevant for the 
description of the ground state of le O and that Vucom provides 
the accurate representation of the relevant nucleon-nucleon in- 
teractions. 



V. CONCLUSIONS 

Through the direct comparison of results for the le O 
ground-state energy obtained using the Vucom interaction 
within the IT-NCSM, IT-CI, and CC approaches, we have es- 
tablished a comprehensive picture of the quality of the dif- 
ferent many-body approaches and the practical relevance of 
formal limitations associated with each one of them. Among 
the points that we have discussed in detail are the possible role 
of CM contaminations, the choice of the single-particle basis, 
and the impact of size-extensivity. 

The analysis of the coupling of intrinsic and CM motions 
using the Lawson prescription shows that the IT-CI and CC 
methods, which are based on a single-particle truncation when 
constructing the relevant model spaces, exhibit a coupling be- 
tween intrinsic and CM motions that in the case of le O af- 
fects the intrinsic energies at the level of about 2-6 MeV or so 
when the HO reference determinant is employed. The same 
analysis also indicates that a small value of the expectation 
value of the CM Hamiltonian (H cm ) alone does not automat- 
ically warrant a decoupling and a spuriousness-free intrinsic 
state. Only the IT-NCSM approach, which is based on ap- 
proximating the complete N max Ml model space of the NCSM 
theory, shows a virtually perfect decoupling, leading to effec- 
tively contamination-free intrinsic eigenstates. On the other 
hand, the relatively small CM contaminations observed in the 
IT-CI and CC calculations for le O do not seem to be detri- 
mental for the quality of the resulting energies, which almost 
perfectly agree with one another when the appropriate levels 
of both theories are employed, namely, IT-CI(4p4h)+MRD 
in the case of IT-CI and CR-CC(2,3) in the case of CC, and 
which are in reasonable agreement with the results of the 
converged IT-NCSM(4p4h)+MRD calculations when we at- 
tempt to extrapolate the IT-CI(4p4h)+MRD and CR-CC(2,3) 
results to the complete basis set limit. The remaining few 
(up to about 10) MeV differences between the crudely extrap- 
olated CR-CC(2,3) and IT-CI(4p4h)+MRD energies and the 
converged IT-NCSM(4p4h)+MRD calculations are consistent 
with the magnitude of the CM contaminations present in the 
CR-CC(2,3) and IT-CI(4p4h)+MRD calculations. 

The comparison of the calculations employing the HO 
single-particle bases with the calculations using HF-optimized 
bases demonstrates that the CC method is largely insensitive 
to the basis set choice, whereas the IT-CI calculations show 
a strong basis set dependence, significantly benefiting from 
the use of an optimized basis. This fundamentally different 
behavior of the CC and IT-CI approaches with regard to the 
choice of the single-particle basis is related to the presence 
of the exp(Ti) component in the exponential wave operator 
of CC theory, as in the Thouless theorem, which makes the 
CC results almost insensitive to the basis set choice, and the 



incomplete treatment of this component by the linear wave op- 
erator of IT-CI. As shown in this work, an effective measure 
of the suitability of the single-particle basis for the truncated 
CI (e.g. IT-CI) calculations can be provided by the T\ or C\ 
excitation amplitudes and the corresponding T\ Diagnostic as 
well as the magnitude of the Davidson extensivity corrections. 
The T\ or C\ values, the values of the T\ Diagnostic, and the 
Davidson extensivity corrections all become large if the ba- 
sis is ill-adapted to the truncated CI calculations of interest. 
We can, therefore, check the suitability of a given basis set 
for the IT-CI and IT-NCSM calculations by monitoring these 
quantities. 

When using the HF-optimized basis in large scale CR- 
CC(2,3) and IT-CI(4p4h)+MRD calculations, we observe an 
excellent agreement of the ground-state energies of le O for all 
values of Ml and e max - Only for the largest e max = 7 space 
and large Ml values do we observe the slightly larger dif- 
ferences between the CR-CC(2,3) and IT-CI(4p4h)+MRD re- 
sults, on the order of 1 to 2 MeV. This excellent agreement be- 
tween the results of the CR-CC(2,3) and IT-CI(4p4h)+MRD 
calculations demonstrates that neither the violation of strict 
size extensivity by the IT-CI schemes nor the violation of 
the variational principle by the truncated CC schemes are 
of major practical concern in nuclear structure calculations, 
since both the IT-CI and CC methodologies are systemati- 
cally improvable through the inclusion of higher-order many- 
body components in the corresponding excitation operators 
(C m components in the case of IT-CI and cluster compo- 
nents T m in the case of CC) and the use of suitable energy 
corrections (the Davidson corrections in the case of IT-CI 
and the corrections due to the effects of higher-order clusters 
in CC). The IT-CI(4p4h)+MRD and the CR-CC(2,3) results 
converge toward somewhat lower binding energies than the 
IT-NCSM(4p4h)+MRD calculations, but, as already pointed 
out, the observed few MeV or a few percent differences be- 
tween the extrapolated e max — > oo IT-CI(4p4h)+MRD and 
CR-CC(2,3) energies on the one hand and the N max — > oo IT- 
NCSM(4p4h)+MRD energies on the other hand are consistent 
with the effects expected from the presence of CM contamina- 
tions in the CI and CC calculations. Based on all of these ob- 
servations, we conclude that all ab initio schemes used in the 
present work — IT-NCSM, IT-CI, and CC — provide pow- 
erful, affordable, and potentially accurate computational tools 
to tackle the nuclear many-body problem. Due to the com- 
plementarity of the IT-NCSM, IT-CI, and CC methods, com- 
parative computational studies using all of these approaches, 
following the analysis presented in this work, may provide a 
comprehensive and precise picture of nuclear structure of p- 
shell nuclei and beyond. 

Finally, it is interesting to compare the results for the 
ground-state energy of le O obtained in the present study with 
the Vucom two-body interaction with the earlier recent CC cal- 
culations using other two-body interactions. The extrapolated 
ground-state energy obtained in the IT-NCSM(4p4h)+MRD 
calculations with Vucom is (—130 ± 1) MeV, i.e. within 
about 2 MeV from the experimental value of —127.6 MeV. 
Although the CR-CC(2,3) results reported in this work and 
obtained with the same potential are not as well converged 
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with the single-particle basis set, the attempt to extrapolate 
the CR-CC(2,3) energies to the e max — > oo limit produced 
the result which implies overbinding, compared to experi- 
ment, on the order of 3-13 MeV. We believe that most of 
this overbinding and the analogous overbinding observed in 
IT-CI(4p4h)+MRD calculations is due to CM contaminations 
that affect the CCSD and the subsequent CR-CC(2,3) calcula- 
tions as well as the IT-CI(4p4h)+MRD calculations, although 
future, more detailed, studies will be required to verify this 
statement. The CR-CCSD(T) study reported in Ref. S us- 
ing the Idaho-A two-body potential [5] produced, after cor- 
recting the result for the effect of the Coulomb interaction, 
an extrapolated ground-state energy of about —109.3 MeV. A 
similar result (approximately, —112 MeV) was obtained with 
the two-body component of the chiral N 3 LO potential jjj, 
which includes the Coulomb interaction directly [48]. Thus, 
with these chiral two-body interactions, the CR-CCSD(T) ap- 
proach underbinds le O by approximately 17 MeV, which was 
concluded to be due to the effect of missing three-body forces 
l48tl . Similar underbinding has also been observed in the CC 
calculations for le O employing the Argonne V18 potential re- 
ported in Refs. 1-441 l45l l4ql . where the explicit inclusion of 
three-nucleon interactions was used to improve the agreement 
with experiment H 0]. A recent CCSD(T) study based 
on the two-body Mowfe interaction presented in Ref. ll52ll 
(CCSD(T) is anoth er non-iterative treatment of T3 clusters 
in CC theory lll25ll . which is less robust and generally less 
accurate than the CR-CC(2,3) approach used in the present 



study |H, m m 0, H 113 H]), has led to an extrap- 
olated ground-state energy of —148.2 MeV, i.e. the binding 
energy overestimated by more than 20 MeV. These compar- 
isons demonstrate that the Kjcom interaction, in contrast to 
many other two-nucleon interactions, requires no or minimal 
contributions from a three-body force in order to provide a 
reasonably accurate description of the ground-state energy of 
16 0. An extension of these investigations to open-shell sys- 
tems and excited states, including other p-shell nuclei, will be 
presented elsewhere. 
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